NACA TN 2693 Q9Q0 



NATIONAL ADVISORY COMMITTEE 


FOR AERONAUTICS 

TECHNICAL NOTE 2693 


A THEORY AND METHOD FOR APPLYING INTERFEROMETRY TO THE 

MEASUREMENT OF CERTAIN TWO-DIMENSIONAL 

GASEOUS DENSITY FIELDS 

By Walton L. Howes and Donald R. Buchele 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 



Washington 

April 1952 


TECH LIBRARY KAFB, NM 





MCA m 2693 


TABLE OF CaNTEirTS 


TECH LIBRARY KAFB, NM 



Page 


STMIARY 1 

mTRODUCTION 1 

PRINCIPLES OF METHOD 3 

Light Path in Flow Field 4 

Optical Path Relation of Interfering Light Waves 7 

Optical Distortion • 11 

Evaluation Equations 12 

Convergence Considerations. 15 

EVALUATION PROCEDURE 16 


NUMERICAL EXAMPLE 17 

SHIFT OF OBJECT PLANE 22 

GENERAL DISCUSSION AND CONCLUSIONS 24 


APPENDIX A - SYMBOLS 26 

APPENDIX B - LIGHT PATH IN FLOW FIELD . . ’ 29 


APPENDIX C - OPTICAL PATH RELATION, 


32 


APPENDIX D - OPTICAL DISTORTION 


36 


references 


37 



KATIOITAL ADVISOEr COMMl'i'i'Kiil FOE AEROKAETIGS 


TECHNICAL NOTE 2693 

A THEORY AND METHOD FOR APPLYING INTERFEROMETRY TO THE MEASUREMENT 
OF CERTAIN TWO-DIMENSIONAL GASEOUS DEHBITI FIEIie 
By Walton L. Howes and Donald E. Buchele 


SUMMARY 

A theory and method are descrlhed for the application of interfer- 
ometry to the measurement of certain two-dimensional gaaeous density 
fields. The theory includes the effects of optical refraction upon the 
ohseryed interference pattern. Exact equations denoting the relative 
density difference corresponding to an observed interference-fringe 
shift and the optical distortion caused by refraction are derived. 
Corresponding approximations in the form of power series expansions are 
developed to permit practical application of the theory. Expressions 
are derived which permit calculation of the values of the power series 
coefficients from experimental data. 

The approximations were applied to available interference data in 
order to determine the density distribution in a boundary layer formed 
by supersonic air flow along a flat plate. Good agreement was obtained 
between the density distribution calculated frcm the interference data 
and that obtained from pressure-probe measurements . 

Tarlous methods for verifying the theory are considered and an eval- 
uation process is outlined. 


INTRODUCTION 

In investigations of hl^ subsonic and supersonic gas flow, where 
compressibility of the gas necessitates the treatment of density as a 
variable, various optical methods (reference l) have been utilized as a 
means of density measurement. This is possible because density variations 
in optical media such as air act in a measurable way on li^t. More- 
over, optical methods of Investigation possess the advantage of peamaltting 
instantaneous recording of the flow without disturbing the flow. Of the 
optical methods, quantitative investigations by the method of interfer- 
ometry have proved most successful. Interferometric investigations of 
flow fields have been conducted primarily with Interferometers of the 
Zehnder-Mach type. The principles of interferometry and the Zehnder-Mach 
interferometer have been described elsewhere (references 2 and 3) . 
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fflie analysis "by means of Interferometry of density distributions 
within optical media is based on the concept of "optical path differ- 
ences" of interfering li^t waTes: the optical-path is defined as the 
integrated product of the physical path traversed by a given wavefront 
and the refractive index along that path. Thus, for convenience, optical 
path differences will be expressed in teims of refractive index rather 
than density. The two quantities are related by the Lorentz -Lorenz law 
(reference 2) . Also, for purposes of clarification, the mathematical 
development is conducted on the basis of the ray theory rather than the 
wave theory of ll^t propagation. This is possible because li^t rays 
are functions of the ll^t waves in that a ray Indicates the direction 
of li^t propagation, which, in nonpolarizing media, is nomal to an 
advancing wavefront. 

In oases such as that considered herein, where the density varies in 
a direction inclined to the direction of propagation of the li^t, the 
mathematical expression of optical path differences is complicated by 
curvature, or refraction, of the light rays. 

Refraction of li^t in gas flow investigations by the method of inter- 
ferometry has been considered by ¥eyl (unpublished), Bershader (refer- 
ence 4), Wachtell (reference 5), Blue (reference 6), and others. An 
optimum wind tunnel span for minimum optical refraction error was deter- 
mined by P. J. ¥eyl in an unpublished analysis. The first two terms 
(constant gradient assumption) of a Maclaurln aeries expansion were effec- 
tively utilized for refractive index. Optical path differences Introduced 
outside the flow field were neglected. This derivation was applied in 
reference 4 to an investigation of supersonic air flow in a channel. More 
recently, Weyl's analysis was expanded by the inclusion of additional 
optical path differences arising outside the flow field and by an expan- 
sion of the Maclaurln series representation of refractive index to several 
terms (see reference 5). A novel method of analysis was used in refer- 
ence 6 to investigate boundary layers in supersonic flow. 

The present analysis attempts to include refraction effects in a 
practical method for evaluating two-dimensional gaseous density fields 
in oases where the greatest space change of density exists essentially 
in one direction. The analysis was carried out at the RACA Lewis labor- 
atory. One coordinate axis of a Cartesian coordinate system is oriented 
parallel to the density gradient. The untaaown refractive index, or 
density, distribution is then expressed as a Maclaurln series with respect 
to that coordinate. The succeeding analysis is essentially an extension 
and modification of the analysis of reference 5. On the basis of the 
succeeding analysis, a procedure for evaluating the data obtained from 
interferograms is described and applied to determine the density distri- 
bution in a boundary layer formed by supersonic air flow along a flat 
plate. Tarlous methods are discussed for verifying the analysis. 
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PEINCIELEB OF 1®3SQD 

A soliematlc diagram, of tke Zehnder-Mach. interferometer apparatus as 
used in the inyestigation of flow-density fields is sh.cwn in figure 1. 

The experimental apparatus consists of a wind tunnel having a rectangular 
cross section hounded on two sides hy plane -parallel glass windows and 
a Zehnder-Mach interferometer positioned with respect to the wind tunnel 
so that li^t (from a monochromatic point source) comprising one of the 
divided collimated ll^t h earns is Incident to the wind tunnel in a direc- 
tion perpendicular to the windows . !The recombined li^t beams are received 
by a camera which is focussed to image sharply seme desired plane, for 
example, plane in figure 1, within the wind tunnel and parallel 

to the windows. When flow is initiated, the resulting nonuniform density 
field is two-dimensional with all density gradients essentially perpen- 
dicular to the direction of the incident li^t. 

With the interferometer in correct adjustment and a condition of no 
flow existing in the wind timnel, an interference pattern consisting of 
either a field of uniform brightness or a series of alternate bri^t and 
dark bands, or fringes, may be observed at the image plane corresponding 
to the plane I-I of figure 1. The phase of interference produced at 
any point of the image plane I-I by the superposition of coherent waves 
is determined by the optical path difference OH) of those pairs of ll^t 
rays which originate at the same point of the li^t source and Intersect 
at the same point in the image plane. 

When flow conditions are established in the wind tunnel, the inter- 
ference fringes observed at the image plane may change their position and 
shape. The fringe shift at each point of the image is measured with 
respect to some reference patteimi, such as the fringe pattern obtained for 
no -flow conditions, or the fringe pattern representing some specific 
region of the flow field, for example, the free-stream region, for which 
the density can be determined. A specific region of the flow field is 
adopted as the frame of reference in the present report. A detailed 
diagram of the i i gbt paths is shown in figure 2. Because the observed 
fringe shifts are attributable to the establishment of flow conditions 
in the wind tunnel, the ll^t path from the source to the initial flow- 
boundary is the same for both no-flow and flow conditions . Therefore, 
the observed fringe shifts may be regarded as a function of only that 
portion of the li^t path which is bounded by the initial flow boundary , 
ai-vi the image plane I-I. Thus figure 2 includes only the region 

of Interest from the initial flow boundary Bj-Bj of the wind tunnel to 
the image plane I-I. For convenience a ray Eg, which circumvents the 
wind tunnel, is shown superimposed upon the wind tunnel. 
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Ll^t Patti in Plow Field 

In o 2 ?der to relate the unlnovm density field to the ohservatlonal 
data secured "by the method of interferometry, the mathematical relation 
between the density field and the observed fringe shifts must be estab- 
lished. This relation may be obtained by treating the evaluation of the 
two-dimensional field as a problem in one dimension whenever the density 
gradient may be treated as a function of only one Cartesian coordinate. 
Consider two Cartesian coordinate ares y and both parallel to the 
density gradient to be investigated and lying in a plane perpendicular to 
the direction of the ll^t path at the initial flow boundary Bj-Bj . 

(All symbols are defined in appendix A.) Let the y origin be located 
at a convenient point along the density gradient and at the initial 
boundary Bj-Bj — in figures 1 and 2. Let the t| origin be in the plane 
Bj-Bj and at that value of y at which any designated li^t ray 
1 = 0,1,2, . . . enters the flow field. Let t] Increase positively in 
the direction of-increas ing density. Let the units of tj and y be 
identical. Then t] and y are related by the transformation 

T[ =|y - y^l (1) 

where y^_ refers to the y-ordinate value at which the li^t ray i 
enters the flow field, and 



( 2 ) 


The y-ooordlnate-serves to locate individual li^t rays with respect to 
the unknown density field whereas use of the T^-coordinate permits the 
mathematical expression of the ray paths in. the simplest form. 


In order to interpret the unknown density distribution in terms of 
the observed Interference pattern, the paths of li^t rays corresponding 
to interfering ll^t waves must be determined. The li^t path is given 
by Fermat’s principle as that path for which the time of transit between 
two points in space attains a stationary value, usually a minimum. Thus 



stationary value 


(3) 


By definition the velocity v of li^t -propagation is related to the 
refractive index ii of the medium by 


= c/v 


(4) 




where c is the velocity of li^t propagation in vacuum. The differen- 
tial ds of physical path length is given by 
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ds 


lOaen eq.uatlon (3) becomea 



1 + 



dz 



dx *! stationary value 


(5) 


( 6 ) 


Let q^uantltles possessing only the subscript i., referring to a given 
li^t ray, represent values of the quantities at the initial boundary 
Bj-Bj of the unknown density field. Let quantities which are starred, 
or subscripted and starred, refer to values of the quantities at the 
final boundary Bjr-Bp of the unknown density field. Then the solution 
of the variation problem (equation (6)), as obtained from Euler's equa- 
tion (appendix B), is 


1 + 



(7) 


which satisfies the boundary conditions 

0 < X 1 X* 
0 < 1) < q* 



= 0 


and possesses a real solution only when p('n) ^ that is, when the 
refractive index is a monotonic increasing function of in which case 
the light ray curves, or refracts, in the xy -plane toward the direction 
of increased refractive index. 


Outside the unknown density field, the refractive index is constant 
within each medium; ienoe, it follows frcm Fermat's principle that the 
li^t paths are strai^t lines. The external path lengths can be deter- 
mined from the geometrical configuration of the experimental setup. 


Equation (?) may be rewritten as the integral equation 
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the solution of which, depends upon the function p. A series solution 
has "been obtained from, equation (?) by the following method (reference 5); 
assume that the refractive index distribution throu^out that region of 
the flew field traversed by any given li^t ray may be represented by the 
Maclaurin series 



(9) 


which converges to the value p, where p is the refractive index of the 
medium, tj is the displacement of the ll^t ray in the wind tunnel with 
respect to the y-oMinate value at which the ray enters the wind tunnel, 


and by definition. 

tv 

1 



“ 7T 

WL, 

^ I=J± 



Let the ray path in the untoiown density field be represented by the 
Maclaurin series 


Oi. 



( 11 ) 


which is assumed to converge . 
first derivative of the series 


Substitution of the series ^ ^ byT)^ and the 
oax^ in the differential equation (?) 


results in an expression in terms of powers of x and containing coef 
ficlents bv and C(j. Following a suitable approximation, the rela- 
tions of the coefficients bv and c<j of like powers of x are 


°0 

°2 

°6 

°0 




1^1 


i V2 


1 T, T, 2 1 , 2.^ 

180 '^1^2 + 40 ^1 ^3 


20-1 
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5040 “1“2 
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( 12 ) 
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vhere Oq = 0 "because t) *= 0 at i = 0. Kius, the path of a given 
light raj i in the field )i is given to a close approximation "bj 

= I ^ 1 ^^ + ^ ^ 1 ^ 2 ^^ + ^ilo ^ 1 ^ 2 ^ + ^ 1 ^ 3 ^ + 

^5040 ^1^2 + 560 ^1^2^3 TT^ ^1 ^4^ X + . . . (13) 

The location of the li^t ray at the final houndary is found "by 

replacing q hy = |yi*-yj_| and x by x*. Details of the deri- 

vation of equation (13) are given in appendix B. 

In addition to the ray path, the angle a, throu^ which the light 
ray is refracted in the unknown density field, is also of interest. By 
means of the substitution 

tan a a (14) 

in equation (?) the angle a may be expressed as 

cos a = — (15) 

r' 


Optical Path Eelation of Interfering Ligjh"t ‘VTaves 

The urnknown density field and the obseived interference pattern are 
related by the optical path difference of interfering light waves. The 
case of a two-dimensional density field contained in a wind tunnel is 
represented by figure 2. In the figure, line E-j_, which intersects points 
Aj_, A2, A3, A4, and P, represents the path of any arbitrary li^t ray 
which traverses the unknown density field. Superimposed on the diagram, 
for convenience, the line E2 (which intersects points B]_, Bg, Bj, B3, 

B4, the focal point F on the lens L-L, and P), represents the path 
of a light ray which circumvents the unknown field and Intersects the 
light ray E]_ at point P in the selected Image plane I-I. Bays E]_ 
and E2 are constructed parallel to the optical axis at the initial 
flow boundary Bj-Bj and are assumed to originate at the same point of 

the Interfercaneter li^t source. According to the latter condition, the 
w ll^t waves associated with these rays are coherent and therefore may 

interfere . 

The parallel construction of li^t rays E-]_ and E2 implies that 
the phase of interference be constant over the image plane I-I, that is, 
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that the interf ercmeter ty initially adjusted for obtaining an "infinite" 
interference fringe. However, it can he easily shown that the succeeding 
analysis applies without modification for any initial fringe adjustment 
and spacing. If an image screen is located arbitrarily at plane I-I, 
then interference of the waves associated with rays E]_ and Eg 
observed at point P. A similar situation applies to all other ray pairs 
emitted from the same point in the ll^t source and intersecting at- 
plane I-I. 


According to the geometrical theory of optical image formation, 
there exists a conjugate object point associated with each image point 
such as P. In the present case it is assumed that the optical layout— 
is such that the object point associated with P is located in the 
vicinity of the flow field. Ihus ray E 2 intersects, or at least 
appears to intersect, ray E]_ somewhere in the vicinity of the flow 
field. 

Suppose that the camera lens system L-L represents a perfect 
optical system in the sense of Gaussian optica . Then, in general, the 
camera lens will produce a perfect point image of a point object when 
the intervening object and image spaces are represented by atmosphere 
under standard conditions; that is, all li^t rays which Intersect the 
object point and traverse the lens will intersect at the same point in 
the image apace. If the lens system is imperfect or if- a medium, such 
as a glass wind-tunnel window, having a refractive index which differs 
from that of standard atmosphere is interposed between the object and 
image, or if both conditions are present then an imperfect image results. 
The resulting image imperfections are called geometrical aberrations . 

For the present, the lens L-L will be regarded as perfect. 

If the wind-tunnel window is neglected and the refractive index of 
the object and image spaces is considered to be that of standard atmos- 
phere, there exists a conjugate object point which is intersected by 
every ll^t ray which intersects E. This object point is indicated by 
the intersection of the backward extension of ray E]_ from point A 3 
with ray Eg. I3ie intersection is at point lying in plane Oj^-0^, 

which represents the apparent Ga^usslan object plane. Similarly, all 
other rays which intersect at the image plane I-I also appear to inter- 
sect at plane Oj^- 0 ^. 

Now consider the interposition of a wind-tunnel window between 
plane O^-Oj^ and the lens L-L. Assuming that the refractive index of 
the surrounding medium remains that of standard atmosphere, there exists 
a real Gaussian object plane displaced from by a distance 

x(t) Independent of angle ag,*- A real object- located in the plane 
Oj-Og will, in the sense of Gaussian optics, appear to be located at 
plane Oj^-0^ when the object is observed through the window. The 
shift x(t) is given by: 



2339 


31 ^ 


NACA ™ 2693 


9 




Jh 


J. 




( 16 ) 


vliere )j,g^ and i-l^ are the refractlTe Indices of atmosphere and the 
window, respectlyely . 


Fezt consider a train of li^t waves generated from point and 

assume that the refractive index throu^out the object and image spaces 
is that of atmosphere. According to the theory of geometrical optics, 
the optical paths traversed by all li^t waves from a given object point 
to the conjugate image point of a perfect optical system are equal. QBius 
the optical lengths defined by points A3, A4, P and ^A^ ®2^ ®3> 

Bg, B^, F, P are equal. Associated with the ll^t waves are wavefronts; 
a wavefront la defined as the locus of points having a constant optical 
distance from Pj^. Since the refractive index is assumed to be constant 
in the object space, the wavefronts about P^^ are spheres . This con- 
dition follows from the definition of optical path lengbh. Consider the 
wavefront centered at ^A and intersecting point A3. It is apparent 

that Hg^P^^ t= ^g_P^B^. Since the optical path lengths from Pj^ to P 
via points A^ and B4 are also equal, it follows that the optical 
path length of ray from A3 to P equals the optical path length 
of ray E2 from Bg to P so that these latter lengths can be 
neglected because they contribute zero OED. !Ehus the entire OED indi- 
cated by an observed fringe shift S at point P is represented by 
the OED of rays Ej_ and Eg in the space bounded by the initial flow 
boundary Bj-Bj and the arc A3B5. Ehe OED of rays Ej^ and Eg, as 
indicated by a fringe shift S^ at point P, is then given by 

Sg^X = ^ ^ ds - BjBg + " ■t) " ^^a ®3% 

J 


where Sg, the fringe shift in wavelengths X of ll^t, is measured 
with respect to the fringe located at P when p. is replaced by Pg. 
By virtue of equations (5) and (7), equation (17) may be rewritten as 




= (h - 






dx - Pj_x* + (AgAg - t) - p^ B3B5 


fl8) 


where, in addition to the substitution of equations (s) and (7), the 
factor P4X* has been subtracted from the first term of equation (17), 
added to the second term, and the terms have been rearranged. Conse- 
quently, the first term of equation (I8) represents the OED due purely 
to the refractive index difference obtained for an unrefracted ray 



10 


NACA TN 2693 


•wtLlch. intersects the hcrimdaiy at that is, at the ordinate 

Talue a 0, or y a The second and third terms of eq^uation (18) 

represent the additional OPD incurred Between the tunnel Boundaries 
Bj-Bj and Bji-Bji as a result of refractlan of ray E-j_, and the fourth 
and fifth terms indicate the OED introduced By the wind-tunnel window 
and the ai 2 ?space Between the window and camera lens L-L^ respectlyely, 
also as the result of refraction. 

The refractive indices of atmosphere and pg of the flow refer- 

ence region^ with respect to which the fringe shifts are to Be measured^ 
are constant so that a change of refractive index within the wind tunnel 
from that of standard atmosphere to that of the flow reference region 
yields a fringe shift Sg given By 

SgX » (Pq - Pa)^:* (19) 

which corresponds to equation (18) with the exception that the refraction 
teims are zero. SuBtraction of equation (l9) from equation (18) gives 

^ px* ^ 

S^X = (p^ " + [17 \ P^ iz - p^x* +^p^(A 2 Ag - t) - pg B 3 B 5 

( 20 ) 

where 3i « Sg - Sg and S^^ represents the oBserved fringe shift in 

wavelengths corresponding to a given li^t ray 1 and measured with 
respect to the fringe pattern representing the flow reference region. 

Let the location of- the real oBJect plane Oj^-Og Be deteimined 
with respect to the internal surface Bj-Bj. of the final wind tunnel 
window By 


X = Zx* 


( 21 ) 


Then it is shown in appendix C By the application of Snell's law of 
optical refraction and consideration of the geometrical configuration 
in figure 2 that the fringe shift Sj^ in equation (20) is given By 
the exact expression 


S^X 


px 

(Pj.-P 5 )x* + p^ dx 

^Jo 


* * 
p^X +PgKX 






a/ 


2 2 *2 

-H^i -^^i ^a 









^ r~T 

A/ Pa +Pj -Pj -P t 

.A/Pa^-H^l^-Pl^ 


( 22 ) 
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•when -fche camera lena eystem is assumed to he perfect. For application 
to practical prohlems, eq.uatlon (22) may he simplified hy mahing suitable 
approximations^ as shown in appendix 0, -with the result 


SiX 


= (Ui - 




dx 


H^x* - Ex*(p* - (23) 


Eq.uatlon (23) is independent of the -wind-tunnel windows, -faheir effect 
upon the ohserved fringe shift thus being negligible. Assuming 

that the series expansions in eq.uatlons (9) and (ll) are convergent, it 
is also shown in appendix 0 "that -bheir substitution in equation (23) 
yields 


X ^ ^ 

M-i - M-5 «= Si „ - 2 

^ v=1ct- 2,4,6. 


LL 




a+ 1 


ov *a _ 
X + K 



,4,6, 


bvCcTV^ 


*CT 


(24) 


where the coefficients Gq-^ may be expressed in terms of the coeffi- 
cients Ocr by eq.'uating the coefficients of like powers of x j,n the 
expression 


E 

rr=2 


'av 


■(t 



(25) 


Optical Distortion 

The refractive index difference - P 5 given by eq.uation (24) 
is associated -with a given li^t ray 1 which enters the flow field at 
■bhe oidLlnate value y^. However, the fringe shift corresponding to 

the ray i may appear to exist at some oi*dlnate value y^ + (Ay)j^ 
other than y^; that is, there may exist a form of optical distortion 
caused by refraction of the ll^t rays in -bhe flow field. For example, 
in figure 2 -bhe li^t rays ®2 appear to originate at the 

apparent object point P^, whereas the refractive index p^, indicated 
by the fringe- shift by virtue of equation (24), exists at the 

y-ordlnate value containing -bhe point A^_ at which ray E-^ enters the 

flow field. Let -bhe difference be-bween -bhe y-ordlnate values contain- 
ing points and A]_ be termed the "optical distortion in the object 

space", denoted by Ay. Then, with respect to rays E]_ and E 2 , the 
distortion may be represented after consideration of figure 2 by 
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(Ay)^ = A2B3 - (A3B3 + B^A^) ( 26 ) 

or 

(Ay)i « I^Ex* - x(t) 4- tj tan a* - |^i)* + y(a^,t)J ( 27 ) 

By use of Snell's lav and equations (15) and (16) it is demonstrated in 
appendix D that equation (27) can he presented in an exact form analogous 
to that of equation (22) as 



When suitable approximations are applied, equation (28) can he rewritten 
to a close approximation as 

(Ay)j^ .= - T)* (29) 

which like equation (23) is Independent of the wind-tunnel windows. By 
virtue of the differential equation (7) of the ray path, the square- 

root term of equation (29) may he replaced hy , where * 1, 

VVil«Tli 

thus yielding 

- < (30) 

Then Ti!f is replaced hy ) 2 :°^ and *, hy the first deriva- 

tlTe of y C£jX*^ (assuming convergence); equation (30) is now expressed 


in its final form 


<7=2 ,4,6, . . . 


Evaluation Equations 

Equations (24) and (3l) relate ohservahle responses of individual 
li^t rays to the refractive index distrlhutlon in the flow field. 
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Eq.u^tlon (24) may te readily transformed Into an expression for the rela- 
tive density change existing at the value at wliloh a given ll^t 

ray 1 enters the flow field hy applying the G-lads tone -Dale approximation 


h *= 1 + kp (32) 

of the Lorentz -Lorenz law, where p Is the gas density corresponding to 
the refractive Indpx [j., and k Is a constant. Ihus, hy substitution, 
equation (24) becomes 


^1~P5 ^ 

P6 kpg 


00 03 






.A. LL 

1 O'— 2 ,4,6,... V =1 ,4,6, . . . 


^v®ov^ 


(33) 




which may be evaluated If the coefficients by are determined. 


Because of refraction of the ll^t rays traversing the flow field, 
■the observed Interference pattern Is, In general, a distorted represen- 
tation of the density distribution •within the field; the distortion of 
the fringe pattern corresponding to any given ll^t ray being given by 
equation (3l) . In order to determine the true density field, the dis- 
tortion effect must be eliminated or at least minimized In same manner. 


In order to evaluate eqiiatlons (33) and (31), the values of the 
coefficients by, Cg-, and c^y must be determined for each ordinate 

value y^^ at -which any given ll^t ray 1 enters the flow field. 
Actually, If -the values of the coefficients by are determined, the 
corresponding values of the coefficients Cq- and Cpy are given 
automatically because -the coefficients Cq- and cov sre functions of 
by according to equations (12) and (25), respectively. The coeffi- 
cients by may be evaluated if they are expressed as functions of the 
derivatives of the observed interference -fringe shifts with respect to 
y. For any given ll^t ray 1, successive differentiations of the 
assumed refractlve-iiadex-dietrlbutlon equation (9) yields equation (lO), 
where derivatives of the coefficients by are zero because the coeffi- 
cients by are constant for any given light ray. Consider equation (24), 
which characterizes any given ll^t ray 1, where and by possess 

fixed values. With respect to the entire manifold of rays comprising the 
li^t beam, p and by (hence c^y) are variables. Thus equation (24) 

may be re-wrltten without subscripts 


1 as 
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CO 00 


V =1 a*= 2 , 4 , 6 , 


TDmC 


v'^av 


CO CO 

v-' 


a + 1 


Vei a^, 4 , 6 , 


’byQ0"y^ (34) 


Differentiation of ectuation (34) v times with respect to y yields 

^ = /2 V ^vCav _*a 


iy 


'* X* dy'' dy'' 


O’ + 1 


yVtssl ct=2^4^6^» 


v=l o-t=2,4,6, 


\°C7V- 


-*cr 


(35) 

Subsequently, it follows from eq.uation (10) that, at any glTen ordinate- 
value yj^, 

\ /d'^S' 


bv 


vl 


/ 


.f-y''/ 


y^yi 


00 00 

r-\ 


b„c 


j^v 1 l—J o 1 

y v=l o»=2,4,6, . . . 


v^ov 3^*0 _ g- 




hh . ^vdjvX 

V =1 ois2,4, 6, . . . 


*o 


(36) 


For simplification, assume that all tenns on the ri^t-hand side of 
eq.uatlon (36) are small cctmpared with the first tem. Then at any 
given ordinate value y=y^ 


i_ A 

Vl 



y=yj_ 


(37) 


In cases where the preceding assumption Is Invalid, the effect of 
the neglected terms can be accounted for by the application of an iter- 
ation process, where "the first approximation of by is given by equa- 
tion (37). 

d'^S 

It- is evident that the derivatives — ^ are determinable from the 
observed Interference pattern. However, i? Is important to note that 

^Vg 

the true values of — hence the desired values of bv; are obtained 

dy'^ 

only when -the distortion Ay is effectively zero -bhrou^out ’the region 
of interest. However, it follows from consideration of eq.uation (3l), 
that, in general, no single object p]ene of focus exists for which 


Ay = 0 for all ll^t-rays. If the ray path series 


E”"' 


■x^ is con- 


vergent, then the first nonzero term - for which o = 2 - may generally 
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te expected to "be the tern of largest ahaolute magnitude . In the case 
of equation (3l) this term may he reduced to zero hy selecting E = 0.5^ 
that Is, hy selecting the object plane of focus to he the midspan plane 
of the -wind tunnel. lEhen eq_uatlon (31) hecomes 


(Ay)^ « ^ ^ t^x**^ (38) 

0 = 4 , 6 , 8 , . . . 

Assuming convergence, the magnitude of equation (38) will usually he 
small enou^ so that the shape of the curve S(y) will not he changed 

d'^S 

appreciably. Thus the values of the derivatives — ^ used to compute 

liy 

the values of the coefficients hv will he essentially unchanged hy 
the residual distortion. If the curve S(y) were changed appreciably 
hy consideration of eq.uatlon (38), then an Iteration process utilizing 
eq.uatlons (37), (36), and (38) could he performed until the desired 
values of hv and the desired curve S(y) were obtained such that by 
satisfied equation (36) and Ay -*• 0 throu^out the Interval 0 1 y < yg 


For E = 0.5, expression (33) hecomes 


Pl"P8 

P5 


1 

kPg 




w 

-^EE 

V=1 a=2,4,6,« 


hyCorv 
a + 1 



09 CO 



v=lais2,4,6, 


^v°ov 



(39) 

Equations (37), (l2), (25), (39), (38), and possibly (36) are suffl 
dent for evaluating an untaiown gaseous density field satisfying the 
restrictions Imposed In the derivation. 


Convergence Considerations 

It Is Impractical to Investigate convergence of the series expan- 
sions and Iteration processes analytically because the exact analytic 

P“Pr 

expressions for Ay and ~ are unknown. The validity of the anal— 

Pg 

ysls can he Indicated hy comparing the results obtained hy the Interfer- 
ence method with the results obtained from other methods of analysis. 

For example, the resultant distribution obtained hy the Interference 
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method can he compared with that given hy the theory for the density 
field or hy some other valid escperimental method. Agreement of the 
series solution with the other methods of solution would Indicate that 
the series are either uniformly convergent or asyniptotic representations 
of the unknown quantities. It is possible that- the series might he 
asymptotic, although divergent, because the Maclaurin series possesses 
the form of an asymptotic series (reference 7) . 

If the series (equations (38) and (39)) are convergent or asymptotic, 
then in computations utilizing these series the remainder in each case 
is less in absolute value than the first series term which is neglected. 
Thus the number of applicable terms of either series will include all 
successive terms for which the coefficients can be calculated with a 
reasonable degree of accuracy up to, but not including, the term of mini- 
mum absolute magnitude. 


IVALUATIOR' PROCEDDEE 


Assume that the interferometer and wind tunnel are in operational 
adjustment and that the interference fringes are obtainable for any 
desired object plane for which the interferometer camera may be focused 
(before flow is initiated in the wind tunnel) by observing at plane I-I 
(in fig. 2) the image formed by li^t reflected from an object inserted 
into the field of view at the desired object plane of focus, usually the 
midspan plane of the wind tunnel, for which K = 0.5. After the camera 
is focussed, the object is removed, flow is initiated, and the interfero- 
grams are obtained by means of photographic processes. 


An interferogram may be evaluated according to any one of several 
schemes depending upon the initial choice of fringe orientation and 
spacing. The result of the interferogram evaluation process is an 
expression for the fringe shift S as a function of y for the partic- 
ular experimental conditions. These values of S(y) can be used in 
the following manner to determine the density field: 


1. From the experimentally determined values of S(y), obtained 

dS 

when E B 0,5, determine the values of successive derivatives gy. 


d^S 


d^'S 


,2’ 


throu^out the interval 0 — y — yg. These values may be 

dy^ dy'' 

determined by any one of several methods j for example, S(y) may be 

fiS\ 

plotted g3c*aphically. Then at any ordinate value y-y^, I5J) 
sents the slope of the curve S(y). The derivatives fg— ] can be 


repre- 


plotted. From the curve values of 


Vdy^ j 


y=7i 

for various values 
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7=j-l can be detertalned as before. Hie process may be repeated to deter- 
mine sucoessiYe deriTatlvea ^ ^ 

dy'' 

Alternatively, S(y) may be represented by a polynomial S = Pv(y)> 


The values of successive derivatives 




at various ordinate 


values jfay^ may then be calculated by taking successive derivatives of 
•with respect to y and then inserting the values j=7±. 

2 . Calculate values of the coefficients by for each li^t ray 
from the values of im 


\<iy Vybsy^ 


and equation (37) . 


3. Determine from the preceding results by(y) the magnitude of 
the expression which is contained in equation (36) but which is neglected 
in equation (37). (The coefficients C(jy are given by equation (25)). 

Calculate the residual distortion given by equation (38), where the coef- 
ficients C{j are functions of the coefficients b^ by virtue of equa- 
tions (12). If the magnitude of the expression neglected in equation (37) 
is of the order of magnitude of the right-hand side of equation (37), or 
if correction of the residual distortion of the curve S(y) results in 


WVy= 


an appreciable modification of the values of the quantities 

■y=yi 

or if both these conditions are true, then the iteration processes men- 
tioned in the section Evaluation Equations may possibly be applied to 
determine more accurate values of (dy)i and b^. 


Pf-Pg 

4. Calculate the relative density difference given by equa- 

P8 

tlon (39) for various ordinate values 7=7^> 

P-P5 

5. Plot against y, where the observed values y^^ are 

replaced by the distortion-corrected values y^^ + (Ay)^^. The resulting 
curve represents the desired density profile. 


HCMEBICAL EXAMPIE 

In order to Illustrate the application of the equations and proce- 
dure and to investigate attendant conditions, consider the evaluation 
of the density distribution along the density gradient within a boundary 
layer due to supersonic air flow along a flat plate. 
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llie particular setup In-veatigated (reference 6) consisted of a flat 
plate of 3. 6-inch width contained in a wind tunnel of 3. 6 -inch span 
hounded hy 0 .5-lnch -thick glass -windows haying a refract iye index of 
1.52. Ihe inteitferometer ll^t source consisted of a high -pres s-ure 
mercury arc in conjunction with an interference-type filter which -trans- 
mitted the green spectral line of wayelength 5461 Angstroms (2.15x10“^ in.). 
Ihe interferometer camera was focussed to image the midspan plane of the 
wind tunnel. 

Supersonic air flow haying a free -stream Mach ninaher of-2.0 at a 
density of 6.77X10“^ slugs per cubic foot was established. Ihe follow- 
ing data were obtained from an interf erogram of -the boundary layer 
2.5 inches downstream of the leading edge of the flat plate: 


Fringe Shift, S 

0 

-.50 

-1.00 

-1.50 

-2.00 

Distance from plate, y. Inches 

0.0429 

0.0198 

0.0165 

0.0142 

0.0120 

S 

-2.50 

-3.00 

-3.50 

-4.00 

-4.50 

y, inches 

0.0100 

0.0083 

0.0066 

0.0046 

0.0028 


Ihese data are presented graphically in figure 3 
yatives 


d% 


^ and 

a5(y) 

resulting curyes — 


The first -two derl- 
The 


were determined -with a tangentiometer . 

^ are plotted in figures 4 and 5^ 


and 


respectiyely. The q^uantltatiye determination of succeeding deriyatlyes 
was of- guest ionable yalue. The residual distortion (indicated by egua- 
tion (38)), the inherent error of the original data oxirye, and the error 
Introduced by -the measuring Instrument would generally- tend to introduce 
large relatlye errors in the measured -values of succeeding derlyatives. 


Ihe measured yalues of ^ and and expression (37) were 

iy dy2 

utilized to calculate- the coefficients b]_ and bg for yarious yalues 
of y to obtain the curyes b^^(y) and t>2(y), which are also plotted 
in figures 4 and 5, respectiyely. 


With respect to the eyaluation eguatlons, if refraction of the 
ll^t rays in the flow field is neglected, then series (9) and (ll) 
reduce to u = bg and T) « 0, respectiyely, and the eyaluation egua- 

tlons (38) and (39) became, respectiyely, 

(Ay)'^ = 0 (42) 


9 


k 


and 
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P 1 -P 6 




Pq i irpgx* 


(43) 


■ffMcli represent first approzimatlons of the residual distortion and 
relatlTe density difference. These equations haye "been utilized hy 
several investigators vith considerable success, especially -when the 
optical refraction vas small. Equations (42) and (43) were used to 
determine a first approximation 'of the density distribution indicated 
by the data of the present example. The resultant distribution is 
shown in figure 6. Then, by means of equations (37), (12), (25), (38), 


dS 


d^S 


and ( 39 ) and the derivatives ^ and — 2 > density distribution 

was calculated to a third approximation. The appropriate approximations 
were determined as follows : 


From equation (37) it follows that 



(44a) 


(44b) 


By virtue of equation (12), equation (38) becomes 

(Ay)i = ^ bj^bgX*^ (45) 

where only terras possessing coefficients byj v < 2,’ hence Cq-, 4, 

d^g 

are Included because only the values of derivatives — v ^ 2, 

determined. On this basis, the expansion of equation (39) is 

Pi"Ps 1 


were 


P 5 kPj 


l^i ^ " i^^l°21 + ^ (^1°41 ''^2°42)^*^ 


However, by virtue of equation (25) 


Cg^x*2 + c^^x*4 = CgX*2 + 


and 
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so tliat 




021 - °2 °41 


°22 “ °42 “ °2 


'Therefore, considering eq.uation (12), the relative density difference 
is given hy 


Pj-P5 _ 

Pg “ kPg 




2 _ 1 . 
i X* “ 12 





) 


(46) 


in th,e present instance. 


The resultant density distribution (third approximation) indicated 
by eq.uatlons (44a), (44b), (45), and (46) is shown in figure 6 for com- 
parison vlth the first approximation. The second approximation is not 
shown in figure 6 because the resulting curve is nearly identical to 
that given by the third approximation ♦ In addition, the third approxi- 
mation of the density profile is depicted in figure 7 for comparison 
vlth the theoretical profile shown in reference 6 and the experimental 
profile calculated in TOferenoe 6 from total pressure measurements. 

The remainders of equations (45) and (46) were obtained by expand- 
ing equations (38) and (39), respectively, to include terms containing 
X*®, with the result 

and 

^ ‘ (li? I’M + jIo (“) 


Considering figure 5, the maximum absolute value of bj was found to be 
about 5 inches"^. Thus the maximum absolute values of the remainders 


were found to be 




K P-P6 


< 0.0001 inch and 

Ay 

P6 


< 0 . 01 . 


In the preceding example, the appropriate evaluation equations were 
developed for the case where derivatives of S of order hl^er than the 
second are neglected. As noted previously, the calculation of values of 
succeeding derivatives would be Impractical in most oases . Thus the 
applicable formulas correspond to polynomials possessing coefficients 


k 


% 


2339 



MCA Hf 2693 


21 


defined at tlie initial "bovindary value x *> 0, y » !Ehe evaluation 

procedure corresponds to that -which vas previously described vi-fch -fche ’ 
exception -fchat the former equations are replaced hy the expressions 
developed in th6 preceding section; -that is, the coefficients are 

given hy equations (44a) and (44h), the residual distortion is given hy 
equation (45), and the density difference at each ordinate value y^ is 
calculated from equation (46). 

In -hhe preceding exazople, the principal basis for verification of 
■bhe validity of "the present analysis -was afforded hy comparison of the 
density field Indicated hy the interferogram -with -that indicated hy 
■theoiy and hy total pressure measurements. Ihe generally good agreement 
of the theoretical ard. probe results vith -fche interferometer result 
■would tend to indicate that the interferogram analysis is -valid and prob- 
ably indicates to a good approximation the true density distribution in 
■fche present instance. 

The maximum indicated density correction due to ligjht refraction 
amounted to a 4-percent increase in the density difference with respect 
to free stream, or about 9 percent of the -total density difference 
across "the boundary layer, as depicted in figure 6. The maximum correc- 
tion occurred nearest -fche surface of the flat plate. This correction 
resulted in considerably better over-all agreemen-fc between the density 
field indicated by the interferogram and that indicated by theory and 
by pressure measurements. 

Ihe Tnn.xiTniTm residual distortion was found from equation (45) to be 
0.0015 inch at a distance y ta 0.015 inch frcaa the surface of "fche flat 
plate. The' maximum absolute remainders due to termination of the series 
were less than 0.0001 inch distortion and 1-percent density change. 

In the present example, the previously mentioned Iteration processes 
were neglected in -fche determination of the coefficients by and -fche dis- 
tortion Ay. The TnaxiTmiTn relative error of -fche coefflcien-bs b'jj 

resulting from neglect of -fchese processes was estimated by means of a 
rou^ calculation "fco be about 10 or 20 percent. 

The present analysis assumes -fchat the flow field is perfectly -fcwo 
dimensional with the density gradient inclined essentially in one direc- 
tion. This condition is, however, only approximated in practice because, 
for example, in -fche Illustrative example considered herein, "end" effects 
(reference 4) se 3 rve to violate the assumed flow conditions. Boundary 
layers are formed along "the s-urfaces of the -wdnd -tunnel wirdows and 
adjacent to the windows and the flat plate. Ixi the example, the effects 
of these de-viatlons from -fche assumed conditions were assumed to be small. 
Ho-wever, preliminary considerations indicate that "end" effect correc- 
tions would result in better agreement of the profile indicated by inter- 
ference measurements wl-fch "that Indicated by theory and by total pressure 
msas-urenients, especially in the region adjacent to the free stream. 
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Because no lens system Is absolutely perfect“lt is to be expected 
that aberrations of the camera lens -«rlll introduce additional optical 
path differences and distortion (reference 8) especially when the magni- 
tude of the optical refraction in the flow field is lar^ and the wind 
tunnel-to-camera lens distance is also great. !Ehis effect is minimized 
by the use of a camera lens which introduces the minimum possible 
geometrloal aberration . In the analysis and in the example, the optical 
path and distortion effects introduced by the camera lens were regarded 
as negligible. 


SHIFT OF OBJECT PLANE 

Additional eYldenoe of the validity of the present analysis in any 
given situation may be obtained by comparing the interference pattern 
observed at an object plane other than that denoted by K « 0.5 with 
that predicted by theory. The shift of object plane introduces a double 
effect upon the observed interference pattern: the phase of interference 

and the observed position of the interference attributed to a given 
li^t ray are shifted because both are ftmctions of K by virtue of 
ec[uations (24) and (3l), respectively. The fringe shift incurred in 
translating the object plane from the midspan plane, corresponding to 
E = 0.5, to another plane corresponding to K = may be obtained from 
the differential of e<iuation (24) j that is, 

<o CO 

dSi = - ^ ^ dE (49) 

v«sl 0 = 2 , 4, 6, ... 

Similarly, the change of distortion is obtained from the differential of 
equation (31); that is. 


<i(Ay)i = E CTCjjX*^ dE (so) 

(7=2, 4, 6, . . . 

The shifts AS^ and A(Ay)^^ due to a shift AE of the object plane 

are determined by integrating eq_uations (49) and (50) between the approp 
riate limits. Thus, for a shift of object plane from E = 0.5 to 
E = Ej_, AS^ and A(Ay)^ are given, respectively, by 

CO CO 

ASi = - X ^ bvca^pc*'^^^ (51) 

v=l 0=2,4, 6, . . . 
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and 

m 

A(Ay)i « (% - |)^ (52) 

0=2,4, 6, . . . 

Agreement of the predicted ouTTe S(y) at E = K-j_ obtained by means of 

eq.iiations (5l) and (52) with the experimental result S(y) would be 
further Justification of the analysis by means of series. It is to be 
noted that in the case where refraction of the li^t rays in the flow 
field is neglected, no change in the Interference pattern is predicted 
when the object plane is shifted. 0?his condition is apparent from 
equations (42) and (43). Also, it is to be noted that the prediction 
process la not reversible ; that is, the interference pattern corres- 
ponding to the object plane for which E = 0.5 cannot be readily deter- 
mined from the interference ' pattern corresponding to the object plane 
for which E = E]_. Ihls situation occurs because true values of bv 
and c^ can be determined to a close approximation only when Ay « 0, 
in which case the plane for which E = 0.5 generally represents the 
optimum choice of object plane. 


A second, set of data S(y) was available for the object plane 
corresponding to E = 1, that is, the flow -boundary plane Bj-Bj nearest 

the li^t source. By means of the values of S(y) for E « 0.5 and 
equations (5l) and (52), the predicted curve S(y) for E = 1 was 
calculated, where, as previously, only terms possessing coefficients 
bv, V <2, hence c^., a < 4, were included because only the values of 


derivatives 
sicxis are 


d''s 


— V <2, were determined. 
3-y 


Thus the appropriate expres- 


AS4 




and 


A(Ay)i = I b^x*^ + I bj^bgx*^ 


with the respective remainders 


and 


Eas »= - + 10 ^1^3) 

®A(Ay) = ^ (1 ^1^2 + I ■^^3^ 


(53) 

(54) 

(55) 

(56) 
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!Ehe masiimm absolute Tallies of the remainders were found to he 




-1 


and 


!e 


A (Ay) 


<10"^ Incdi. 


The resulting curye S(y) for S e i predicted by the present 
theory is depicted in figure 8 for comparison with the curve obtained 
from, experimental measurem.ents . Althou^ the agreement of the two 
curves is not perfect, it is still considerably better than that obtained 
■when refraction was neglected, as is evident from figure 8. Apparently 
the disagreement between the experimental result and that predicted by 
the present theory is caused principally by the fact that the errors of 
the calculated values of the coefficients by were relatively large when 

the approximation (equation (37)) of equation (36) was used alone. Thus, 
because the coefficients b^, v > 0, comprise the major factors of equa- 
tions (51) and (52), the relative errors of the quantities ASj^ and 
A(Ayj^) were also large. In the case of the evaluation equation (39), 
the effect of errors of the calculated values of the coefficients by 
appears only in correction terms which are small relative to the primary 
term. Of lesser importance is the fact that the values of the derlva- 

d'^S 

tlves — ^ were determined from a curve S(y) which possessed residual 
dy 

distortion. Thus, the measured values of the derivatives would have been 
in error wherever the slope of the experimental cuire S(y) deviated to 
any great extent from the slope of the desired distort ion -free curve. 


It is concluded that the improved agreement of the predicted curve 
S(y) for E = 1 (when refraction of the light rays in the flow field 
is considered) with the experimental ciarve is further evidence of the 
validity of the series analysis, at least for the present example. 


GEEEEAL DISCUSSION AND CONCLUSIONS 

The present analysis represents an extension and modification of 
Wachtell's analysis with the inclusion of a procedure for analyzing 
certain two-dimensional gaseous density fields by means of optlcab- 
interf erence . 

The' equation of the li^t path in the untnown density field was 
developed as a power series, the coefficients of which are functions of 
the coefficients of a Maclaurin series representation of refractive 
index. An exact optical path difference equation was derived, which, * 

when re-expressed in the foim of an approximation, permitted the deter- 
mination of the unknown density distribution from obseirved interference 
fringe shifts . A form of optical distortion resulting from refraction 
of the ll^t rays in the unknown density field was considered. The 


2339 



2339 


5R 


NAG A TW 2693 


25 


distortion was escpressed "by an exact eq^uation and a power series. The 
power series was utilized to determine an object plane of focus for 
which the distortion is jninlmlzed and to calculate the residual distor- 
tion. Differentiation of the optical-path -difference equation permitted 
direct calculation of values of the fiist few refractive index series 
coefficients . 

In addition, a procedure for applying the present theory was dis- 
cussed and applied to determine the density distribution in a boundary 
layer formed by supersonic air flow along a flat plate. The resultant 
distribution agreed reasonably well with that given by theory and by 
pressure -probe measurements . Further confirmation of the theory was 
obtained by predicting relatively successfully the Interference pattern 
produced for an object plane other than the optimum plane. 

It is concluded that the present analysis will permit increased 
accuracy in the determination of certain two-dimensional gaseous density 
distributions by the method of optical Interference. 


Lewis Fli^t Propulsion Laboratory 

National Advisory Oonmittee for Aeronautics 
Cleveland, Ohio, January. 3, 1952 
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APPEUDIX A 
SYMBOLS 

The folloving symbols are used in this report: 

Ai,A 2 ^A 3 ,A^ points intersected by ll^t ray which traverses flow 

field 


Bv 

bv 

c 

Off 

F 

I-I 

K 

k 

L-L 

°A"°A 

Or-Or 

OPD 

P 

Pa 

p(y) 

R 


final flow boundary traversed by li^t ray 
initial flow boundary traversed by li^t ray 
binomial coefficients 

points Intersected by light ray which clrcmvents flow 

refractive index coefficient 

velocity of light in vacuum 

light path coefficient 

focal polnt^f camera lens; function 

image plane 

fraction of wind tunnel span 

constant (0.1166 ft^/slug for air when 5461A) 

camera lens 

apparent object plane 

real object plane 

optical path difference 

image point 

apparent object point 

polynomial 

remainder 
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S 

s 

t 

V 

X 

x(t) 

x^y 

y(oJ,t) 

yhH,4 

a 

A 

\ 

v,a 

p 

Subscripts: 

a 

1 


light rays 

interference-fringe shift (wavelengths of light) 
distance measured along path of ll^t ray 
wind-tunnel window thickness; time 
velocity of light in medium 

distance from final boundary of flow field to real 
Gaussian object plane 

distance from real Gaussian object plane to apparent 
Gaussian object plane 

coordinates 

A 

displacement of ll^t ray in final wind-tunnel window 

denotes value for given light ray at initial flow- 
boundary traversed by ray 

denotes value for given light ray at final flow- 
boundary traversed by ray 

angle measured with respect to x coordinate 

increment (Ay = distortion) 

coordinate parallel to y coordinate 

wavelength of ll^t 

Integers 

refractive index 
density 

atmosphere 

Integer, refers to a given ll^t ray 
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w vlnd-tumel window 

6 flow reference region (free stream) 

integers 

Superscripts: 

* refers ter value at final flow-boundary traversed by 

li^t ray 
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APPiEKDIX B 


LIGHI PATH Bf FLOW FIELD 
The li^t path in the flow field is obtained as the solution of 


the variation problem 

^ " stationary 

where m- = M-(y) • Equation (6) is of the form 

/it ta 

the solution of lAiich is given by the Euler equation 

Sf 

F - Ti' e constant 

where t]’ * solution of equation (6) is 


value 


( 6 ) 


(Bl) 


(B2) 


where |i = at y = y^ 
Maclatirin series 


. Let the light path be expressed by the 


( 7 ) 


CT=0 


( 11 ) 


where the coefficients C(j are given by 




(B3) 


It follows that 


dq _ 
dx “ 



acjx 


. 0-1 


(B4) 
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Hie- as Burned refractive Index distribution 



v=0 


( 9 ) 


expression (B4), arid finally series (ll) are substituted in the- dif- 
ferential equation (?) to obtain the expression 



The ri^t-hand side of equation (BS) may be replaced by its binomial 
expansion. However, only the first two terms of the three-term expan- 
sion need be retained because the coefficient bQ represents the 

refractive index at ti = 0. Since p. is very nearly unity over a 

wide range of tenperatures and pressures, the third term of the binomial 
expansion is negligible. Therefore, equation (BS) becomes 



(B6) 


where bQ * 1. When coefficients of like powers of x are equated, 
the coefficients c^j may be expressed as functions of the refractive 
index coefficients b^. The first few coefficients c^. are 


Cq = Cl = C3 


CCT 5S * • • S O 


20-1 


= 0 


=2 = I ^1 


°4 = ^ 1^2 


( 12 ) 


=6 = So + « ”1^ ”3 


Co = 


bnb 


8 ” 5040 560 ‘"I ^2^3 H2 ‘^l 
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Replacing the coefficients c^y in equation (ll) by the corresponding 

expressions in terms of coefficients b„ results in the following 
expression for the ray path: 

1] = i bn + — bnboX^ + ( — i- bnbo^ + — bn ^ b 3 ^ X® + 

I 2 1 12 ^ ^ V180 40 1 3y 

(solo ^1^2^ + 560 ^1^ ^2^3 + 3J2 ^4:) x® + . . . (l3) 
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APPENDIX C 


OPTICAL PATH RELATION 
Given the fringe shift e3q>ression 

X* 


S^X = (Hi-iig)x* \x^ 6x - -t li^ (AgAj-t) - lia ®3% (^O) 


the last two terms may be expressed hy consideration of figure 2 as 


IV O^-t) = ivt - ij 


( 01 ) 


and 


Ma B 3 B 5 - - IV [k'O' - f 


respectively, where 


X = 


( 21 ) 


The factors and a* in equations (Cl) and (C2), respectively, 

may he expressed in terms of a* hy virture of Snell’s law of optical 
refraction; that is. 


^L* sin a*- s sin ot* = |ig^ sin a* 
and x(t) may he replaced hy 

-M = 

Then equations (Cl) and (C2) become, respectively, 

1 


(AgAj-t) = ti^t 




(C3) 


(16) 


(04) 




2 i 
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axid 


^^a B3B5 = - [ia ^ 


1 - 




(C5) 


sin^ a* 


G?heii, when (l-cos^ a*) Is substituted for sin^ a* and cos a* is 
replaced by expression (is), that is. 


cos a* = 




equations (C4) (cs) may be rewritten in the form 

tv 


tv (AgAs-t) = 


+ t^i^ - 


- 1 


(C6) 


and 


t% B3®5 “ - tia + 


-t'l 

tv / 


1 - 


A /„ 2 , ,, 2 ^2 

'V'^a ^^i " *^i 


(C7) 


respectively. Substituting equations (C6) and (C7) in equation (20) 
and collecting terms yield 


S^X 


(tJ’i-tis) 35* + 


_1 

t^i 



dx - 



t 


tv 


tv 


-4 


.. 2 . 2 

hV + t^i 


Vtv 


2 + ij .^2 _ ^*2 




( 22 ) 


The ri^t-hand side of expression (22') consists of a sequence of terms 
comprising differences of quantities of nearly equal value. Those terms 
which are included as roots may be ej^ianded by means of the binomial 
theorem. Thus 
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1 1 

■where the hinomlal coefficients By are Bq = 1, B^_ » Bg = - - q , 

. . . The series (C8) converge hecause ^^2 >> ^*2 _ ^^^2^ 

Moreover, because for v >i the terms of expansions (C8) conprise pro- 
ducts of- differences of quantities of nearly equal value, such terms 
may be regarded as negligible. Hence- expression (22) may be reduced to 

X* 

S^X = (ni-ti5)x*-+ ^ C - li^x* - Kx*^ja|2_^^2j 

where li^* 1 except when differences of q'uantities of nearly equal value 
are concerned. Furthermore, because 

_ ^^2 « -ii^) (CIO) 

so that approximation (C9) may be re'written as 

X* 

Si'- =• + ir / 1“^ '3^ - -l^l) (S3) 


Substituting expansion (9) in (23) and replacing i) by series (ll), 
both expansions being assumed con'vergent, yield the following expression 
for approximation (23); 
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where the third term of the hlnomlal e 2 g)anslon of j_g negllglhle. 

Equation (Cll) may he rewritten in the form 



•where the coefficients are e3cpressed in terms of the coefficients 

Cq hy equating the coefficients of like powers of x in 'the 
expression 



Integrating eq-uation (C12) and solving for - pg result in 


Pi - Pg = Si 3^-2, 






a=2/¥76, . 


V sal 0=2 




(24) 


■where the multiplicative constant Pj_ in the last sum of equation (C12) 
is regarded as unity. 
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APPENDIX D 


OPTICAL DISTORTION 
From equation (27), that- is, 

(Ay)i = [kx* - x(t) + tj tan a* “ jjli + y(ov^'fc)J (27) 

it follows hy virtue of Snell's law (equation (C3)) and equation (is) 
that 

Bln <4 - 1 

<=“s “a ■= -j^ (Dl) 


tan a* =/ 

SL 


*2 2 
^^i - ^^i 


Also, 


y(o*,t) *= t tan a* = -^ t 


sin a* 




sin2 <4 


Then, because of equations (D1), 


y(o>?^;'t) = t/ 


- jx^2 


(D2) 


(D3) 


When x(t) is replaced by expression (l6) and equations (D1) and (D3) 
are utilized, equation (27) becomes 


2339 
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As miming that * iig^» 1 , except where differences of qiiantlties 

of nearly equal value coirprlse a congilete term, equation (28) reduces to 

(Ay)i ^ A/^if - ^29) 

Rather than expand the square root term of approximation (29) by the 
binomial expansion, "which converges very slo"wly because ti* « Hj_, the 

term may be replaced by , "where 1, which follows from 

equation (?) so that 

(^)i - (i)„_, - ’ll 

^ ‘'n=Tii 

Finally, because of series (ll) and its derivative, series (B4), 
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Figure 4. - Varlanlott of Interfexence-frings aJLcpa anil coefficient b]_ vith distance frcn surface 
of flat plate. Traction of vlnd tunnel span, K, 0..5. 
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Flifure 7. 


- Coa^arlson of density profiles indie 







